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ABSTRACT 

If the mean properties of clusters of galaxies are well described by the entropy-driven 
model, the distortion induced by the cluster population on the blackbody spectrum of 
the Cosmic Microwave Background radiation is proportional to the total amount of intr- 
acluster gas while temperature anisotropies are dominated by the contribution of lO^^M© 
clusters. This result depends marginally on cluster parameters and it can be used to 
estimate the number density of clusters with enough hot gas to produce a detectable 
Sunyaev-Zel'dovich effect. Comparing different cosmological models, the relation de- 
pends mainly on the density parameter Q^n- If the number density of clusters could be 
estimated by a different method, then this dependence could be used to constrain fi^- 

Subject headings: Cosmic Microwave Background. Cosmology: theory. Galaxies: clus- 
ters: general. 
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1. Introduction. 

Clusters of galaxies can be detected in the X-ray 
band due to the emission of the intrachister (IC) gas. 
This gas is hot enough to change the brightness of 
the Cosmic Microwave Background (CMB) photons 
through inverse Compton scattering (Zel'dovich & 
Sunyaev 1969). The Sunyaev-Zel'dovich (SZ) effect 
caused by individual clusters has been measured for 
tens of clusters (see Birkinshaw 1999 for a review). 
Eventually, the PLANCK satellite will produce an 
all-sky catalogue likely to contain thousands of SZ 
sources (da Silva et al. 1999). Together with the ef- 
fect of single sources, the overall cluster population 
induces distortions and temperature anisotropies on 
the CMB. The Sunyaev-Zel'dovich effect has a well 
known frequency dependence that helps to distinguish 
it from other foregrounds and methods have already 
been proposed to separate and measure this contribu- 
tion (Hobson et al. 1998). 

Extensive theoretical work has been devoted to an- 
alyze the effect of clusters on the CMB radiation (Cole 
& Kaiser 1988, Bartlett & Silk 1994, Colafrancesco et 
al. 1994, Atrio-Barandela & Miicket 1999 -hereafter 
paper I-, Komatsu & Kitayama. 1999). Barbosa et 
al. (1996) noticed that the value of the mean Comp- 
tonization parameter, that measures the amplitude 
of the blackbody distortion, depends on which are 
the less massive clusters that produce a significant 
effect. On the other hand, Cole & Kaiser (1988) 
remarked that temperature anisotropies were domi- 
nated by massive clusters at moderate redshifts. Since 
clusters of a given mass contribute differently to tem- 
perature anisotropies and distortions, in this letter we 
show how the number density of clusters with enough 
hot gas to produce a detectable SZ effect can be esti- 
mated. This number is marginally dependent on the 
cluster model and it varies by at most a factor of four 
in different cosmologies. 

2. Scaling relations of clusters. 

Following paper 1, we assume that clusters are 
spherical in shape, with typical size the virial ra- 
dius r^, virial mass M and have virialized at red- 
shift z. We also assume that the IC gas is isother- 
mal, distributed smoothly and with a radial distribu- 
tion well fitted by the spherical isothermal (3 model 
ne(r) = nc[l + [r/rcf]-^'^/'^ (Jones & Forman 1984), 
where the central electron density and the 

core radius of the cluster, respectively. For conve- 



nience, we shall assume that (3 = 2/3. Similar results 
to those presented here are obtained for any value 
of (3 within the observed range (0.5 < /3 < 0.7, see 
Markevitch et al. 1997). Further, we assume that the 
dynamical evolution of the IC gas is well described by 
the entropy-driven model by Bower (1997) since this 
model seems to provide an adequate description of 
clusters with redshifts 2 < 0.4 (Mushotzky & Scharf 
1997). 

To compute the spectral distortion and tempera- 
ture anisotropies induced by hot IC gas we need to 
translate the properties of a sample of clusters at low 
redshifts into their equivalents at earlier epochs. For 
the spherical collapse model, the virial radius scales 
with mass and redshift as = ryo{M /IQ'^^Mq)'^/'^{1+ 
z)~^ . If the electron temperature is proportional 
to the velocity dispersion of the dark matter then 
T = Tg(M/lOi5Aif0)2/3(l + z). In the entropy-driven 
model, the electron central density scales as ric = 
nco{T/Tgf/^{l + z)-3<^/2. Finally the core radius 
scales as Tc = reo{M/l(}^^MQ)-^/^{l + zf-'^+^''y'^. In 
this scaling relations, the current 

average core and virial radius, the central gas tem- 
perature and electron density of a IO^^Mq cluster, 
respectively, while e parametrizes the rate of core en- 
tropy evolution. Mushotzky and Scharf (1997) found 
e = 0±0.9. 

3. Distortions and Temperature anisotropies. 

The effect of the IC gas is both to distort the black- 
body spectrum and to induce temperature anisotropies 
on the CMB. A measure of such distortion is the mean 
Comptonization parameter defined as 

y = gikf/KsTo) J ^dM^dzKy^^. (1) 

In this expression g{x) = a;coth(a::/2) — 4 gives the de- 
pendence of the SZ effect with frequency /; To is the 
CMB mean temperature, dn/dM is the cluster num- 
ber density per unit of mass, r/o = [ksiyT /'meC^]f'cTgnc 
with physical constants having their usual meaning 
and K — {ry/2dA)'^ gives the probabihty that a par- 
ticular line of sight crosses a cluster, with dA the angu- 
lar distance to the cluster. If we introduce p = ry/vc 
then (f> — Ap^'^{p — tan~-^p) is the averaged line of 
sight through a cluster. 

The effect of all clusters on an angular scale I can 
be obtained by adding in quadrature the contribution 
of each single cluster. If we assume that clusters are 
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Poisson distributed on the sky the radiation power 
spectrum is (Cole & Kaiser 1988) 

Pil) = I ^dM^dzig{x)y,f\m': (2) 

where 4>(l) is the Fourier transform of the angular pro- 
file of the IC gas. This expression represents the con- 
tribution of the foreground cluster population to the 
power spectrum of CMB temperature anisotropics. 
The variance of the temperature field is given by 

AT 1 /■ 

< (^)' ^ y idiw{i)p{i) (3) 

In this equation, W{1) represents the window func- 
tion of the experiment. This temperature anisotropy, 
together with other secondary contributions, adds in 
quadrature with the "intrinsic" (cosmological) signal 
to give the pattern of temperature anisotropics on the 
sky (see paper I for details). 

To compute the number density of clusters as 
a function of mass and redshift we use the Press- 
Schechter approximation, known to describe accu- 
rately the abundance of clusters in numerical simu- 
lations at different redshifts 

dM V TT^M' a dM ■ \ I 

In this expression, v = (Syb/as) is the peak-height 
threshold, a is the rms mass fluctuation within a top- 
hat filter measured in units of erg , the rms mass fluctu- 
ation on a sphere of 8/i~^Mpc today, is the thresh- 
old overdensity of the spherical collapse model and b 
is the bias factor. 

4. The cluster number density. 

We integrate numerically eqs. (|l|) and (||) for three 
cosmological models with the same baryon fraction 
Qb = 0.05: two flat models, (1) standard CDM, 
with Qcdm = 0.95, (2) ACDM with f^cdm = 0.25 and 
f^A = 0.7 as indicated by the recent measurements 
of high redshift supernovae (Perlmutter et al. 1999) 
and (3) an open model, OCDM, with Qcdm = 0.25 
and no cosmological constant, consistent with clus- 
ter dynamics and evolution (Bahcall 1999). As aver- 
age parameters, for a IO^^Mq cluster we took: r^o = 
1.3/i-iMpc, Tg = lO^K, rco = O.Uh-^Mpc and Uco = 
1.19 • 10~'^/i^/^cm~^. We considered different gas evo- 
lutionary histories, by taking e — —1,0,1. The matter 



power spectrum was normalized to cr(8/i~^Mpc) — 0.7 
(Einasto et al. 1999b) and we considered the follow- 
ing bias factors: b — 1.05, 1.3. The first corresponds 
to the bias of QDOT galaxies and the second is an av- 
erage value obtained from several catalogues (Einasto 
et al 1999a). Also, we considered 6y = 1.5, 1.7 (Tozzi 
& Governato, 1998, Governato et al. 1999). The sim- 
plicity of the Press-Schechter approach lies on the fact 
that the number density depends only on the ratio 
v = Syb/as- With the values quoted above, v varies 
in the range 2 to 3 and we shall quote our results for 
both limits. 

To compute the effect of the cluster population on 
the CMB radiation it is necessary to determine first if 
groups of galaxies of about lO^^M© have enough hot 
gas to produce a significant contribution. The scaling 
relations of Sec. [2] have been found to be accurate for 
clusters above 1O^^M0. However, they could very well 
be extended below this mass scale. In this respect, it 
is necessary to specify the integration range of eqs. (P 
and (^). We tried three different lower mass limits: 
Mrmn = 0.1,0.5 and 1 (in units of IO^^Mq). The 
main conclusion was that when changing the mass 
limit one order of magnitude y could vary a factor 3 
to 10, depending on the model, but the temperature 
anisotropy AT/To varied by less than 30%. In this 
calculations and in the rest of the paper, temperature 
anisotropics were computed using the window func- 
tion of the SuZie experiment (Church et al. 1997). 
For the models considered above, the anisotropics 
in the Rayleigh- Jeans regime ranged from 14/^K for 
sCDM with = 2 to l^K for ACDM with ly = S. 

Integration of eq (|^) shows that PP{1) reaches a 
maximum around I = 1000 — 3000. At those scales 
Komatsu & Kitayama (1999) have concluded that the 
increase in temperature anisotropy induced by cluster 
correlations is negligible. To illustrate the different 
behavior of temperature anisotropics and distortions 
with the lower mass integration limit we shall con- 
sider the behavior of the radiation power spectrum at 
I = 1000. This scale is close to the maximum of the SZ 
effect and is a good compromise between the smallest 
angular scale that will be resolved by PLANCK and 
the maximum of the SuZie window function. In Fig. 
1 we show dy/dM, dP{l = lGGG)/dA/ for different cos- 
mological models. Notice that the integrand of y in- 
creases at the low mass end, but that of the radiation 
power spectrum reaches a maximum around 10^'^Mq. 
Its exact location is weakly dependent on model pa- 
rameters. A qualitative explanation of Fig. 1 can be 



3 



given: since tr ~ 1 on a mass scale of « 5 x 10^ M©, in 
the upper mass limit dy/dM and dP/dM are damped 
by the exponential factor in the Press-Schechter for- 
mula; in the low mass limit, the exponential fac- 
tor is close to unity and dy/dM ~ Af(i+")/6 but 
dP/dM ~ M^/^^+^Z^ where n is the slope of the 
matter power spectrum at lO^^M©, roughly n ~ —2 
(Einasto et al. 1999a). 

To conclude, the gas in groups contributes signif- 
icantly to Comptonization but little to temperature 
anisotropics. This fact can be used to estimate the 
number density of clusters that produce a significant 
SZ effect. For this purpose, in paper I we introduced 
the parameter ij = ^-^-j^. It represents a relation 
similar to the approximate solution of the full kinetic 
equation for the change of photon distribution due 
to inverse Compton scattering for a single cluster: 
at/To = g{hf /kTo)yo- As an example, in Fig. 2a 
we plot 7] with respect to lower mass integration limit 
for three cosmological models with e = 0: solid line 
(sCDM), dot-dashed Une (OCDM) and dashed Une 
(ACDM). Thick lines correspond to v — 3 and thin 
lines V = 2. Similar behavior can be observed for 
other gas evolution histories and model parameters. 

By inspection of eqs. (|l|) and (||), the Comptoniza- 
tion parameter and the power spectrum scale linearly 
with the number density of clusters. Therefore, one 

— 1/2 

expects that -q will scale as n^j . In fig. 2b we 

1/2 

plot r/n^i for the same models as in Fig. 2a. No- 

1/2 

tice that the value of rjn^i is almost constant and 
nearly independent on the threshold i/. The actual 
value varies by a factor of two when the total mat- 
ter content changes from fi^ = 1 to 0.3 but does not 
depend on the geometry of the cosmological model, 
i.e. on the value of A. In Fig. 3 we plot as a 
function of gas evolution history, assuming that the 
number density of all the objects contributing to the 
SZ effect is 2.5 x IQ-^h^Mpc^^. In each panel, thick 
lines correspond to a threshold v = S and thin lines 
to = 2. Solid lines correspond to a lower mass 
integration limit of lO^'^M© and dot-dashed lines to 
lO^^M©. Notice that r] varies by less than 30% for dif- 
ferent gas evolutionary histories. Let as remark that, 
if rim = ils + ilcdm is kuown, T] permits to obtain 
rici- Vice versa, if rid could be determined by other 
means, rim could be determined. If no assumption 
is made about the cosmology, rid can be determined 
within a factor of four. 



5. Discussion. 

The results of previous section indicates that rj is 
mostly sensitive to the number density of clusters, and 
for a given cosmological model is weakly dependent on 
the exact modeling of the cluster population. There- 
fore, a measurement of both the mean distortion of 
the CMB radiation and the temperature anisotropics 
induced by clusters of galaxies will permit to estimate 
the number density of objects that have enough hot 
gas to produce a detectable effect on the CMB. This 
conclusion is not limited by the validity of the scal- 
ing relations of Sec. [2]. One must expect a similar 
behavior for a different cluster model, even though 
the level of temperature anisotropics and distortions 
will be different. If the hot IC gas was uniformly dis- 
tributed in the Universe, it would still produce black- 
body distortions but will not give rise to temperature 
anisotropics, and rj would be zero. Only when the 
gas is clumped in clusters, two different lines of sight, 
seeing different column densities, will have different 
temperatures. Therefore, y scales with the number of 
clusters along the line of sight while for Poisson dis- 
tributed clusters AT /To scales as the square root and 

— 1/2 

rj oc , as shown. Therefore, independently of 

the exact model used to describe the average cluster 
properties, 77 measures how the hot gas is distributed 
in the Universe. Finally, if the cluster number density 
is known, 77 could be used as an estimator of the total 
matter density flm- 

To summarize, clusters of IO^^Mq mass dominate 
the contribution to temperature anisotropics. There- 
fore, even if the adequacy of the scaling relations of 
Sec. [2] below that mass scale has not been estab- 
lished, temperature anisotropics can be reliably cal- 
culated using eqs. (||) and (||). If the entropy-driven 
model proves not to be adequate to describe groups 
of galaxies, one would have to use hydrodynamical 
simulations such as those of da Silva et al. (1999) 
to compute y but 77 would still provide an estimate of 
the number density of SZ clusters. Comparison of this 
density with that of X-ray clusters will undoubtably 
help our understanding of cluster formation and evo- 
lution. 
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Fig. 1. — Integrand of the mean Comptonization parameter (curves on the left side) and the radiation power 

spectrum (on the right) for different models. Solid lines correspond to standard CDM. dashed lines to ACDM and 
dot-dashed line to OCDM. The scale in the y-axis is arbitrary; the data has been rescaled for convenience. The 
amplitudes for dy/dM at lO^^Mg and dP/dM at the maximum are: 3 x 10"^, 1.1 x 10"^^ for sCDM, 1.5 x 10""^, 
1.4 X 10-1^ for OCDM and 1.3 x 10"^, 1 x 10"^^ for ACDM. All these quantities are in units of {M/W^Uq)-^ 
and we took g{x) = 1. 
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Fig. 2. — (a) Variation of tj (three upper curves) with the lower mass hmit of integration Mmin for three different 
cosmological models: SCDM (solid line), OCDM (dot-dashed line) and ACDM (dashed line). Thick lines correspond 

1 /2 

to a threshold v = 3 and thin lines to v = 2. (b) Variation of rjnJi a function of Mmin for the same cosmological 
models as in (a). 
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Fig. 3. — The value of ry rescaled to a cluster density of 2.5 x 10~'*ft-^Mpc~^ clusters, as a function of the gas 
evolutionary history. From left to right: sCDM, OCDM and ACDM. Thick lines correspond to = 3 and thin lines 
to z/ = 2; solid lines to Mmin = 1 and dot-dashed lines to Mmin = 0.1 (in units of lO^^M©). 
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